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Abstract. Sudden turn-on of a matter-wave source leads to characteristic oscillations 
of the density profile which are the hallmark feature of diffraction in time. The 
apodization of matter waves relies on the use of smooth aperture functions which 
suppress such oscillations. The analytical dynamics of non-interacting bosons 
arising with different aperture functions are discussed systematically for switching- 
on processes, and for single and many-pulse formation procedures. The possibility 
and time scale of a revival of the diffraction-in-time pattern is also analysed. Similar 
modulations in time of the pulsed output coupling in atom-lasers are responsible for 
the dynamical evolution and characteristics of the beam profile. For multiple pulses, 
different phase schemes and regimes are described and compared. Strongly overlapping 
pulses lead to a saturated, constant beam profile in time and space, up to the revival 
phenomenon. 



PACS numbers: 03.75.Pp, 03.75.-b, 03.75.Be 
Coherent matter-wave pulse formation and dynamics have been studied traditionally 
in the context of scattering and interferometric experiments with mechanical shutters. 
Suddenly switching-on and off matter-wave sources leads to an oscillatory pattern in 
the particle density which was discovered by one of the authors in 1952 and dubbed 
diffraction in time [1, 2] because of the analogy with the diffraction of a light beam from 
a semi-infinite plane. The (sudden) "Moshinsky shutter" describes the evolution of a 
truncated plane wave suddenly released and admits an analytical solution which remains 
a basic reference for analysing more complex and realistic cases [3, 4, 5, 6, 7]. A wide class 
of experiments have reported diffraction in time in neutron [8] and atom optics [9], and 
with electrons as well [10]. It was also observed in a Bose-Einstein condensate bouncing 
off from a vibrating mirror, proving that the effect survives even in the presence of a 
mean-field interaction [11]. A more recent motivation for studying matter- wave pulses 
is the development of atom lasers: intense, coherent and directed matter-wave beams 
extracted from a Bose-Einstein condensate. One of the prototypes uses two-photon 
Raman excitation to create "output coupled" atomic pulses with well defined momentum 
which overlap and form a quasi-continuous beam [12, 13] (For other approaches see 
e.g. [14] or [15]). Clearly, the features of the beam profile are of outmost importance 



Time modulation of atom sources 



2 



for applications [16, 17, 18], and will depend on the pulse shape, duration, emission 
frequency, and initial phase. While many studies on matter-wave pulse formation and 
dynamics deal with single or double pulses, more research is needed to understand and 
control multiple overlapping pulses relevant for current atom-laser experiments. We shall 
undertake this objective within the analytical framework of the elementary Moshinsky 
shutter, here at the level of the Schrodinger equation for independent bosonic atoms. 
Non-linear effects will be numerically considered elsewhere. This formulation may indeed 
be adapted also to other important features of pulse formation in atom lasers, which 
are different from the mechanical shutter scenario. Instead of sudden switching we 
shall consider smooth aperture functions, i.e., "apodization" , a technique well-known 
in Fourier optics to avoid the diffraction effect at the expense of broadening the energy 
distribution [19]. Manifold applications of this technique have also been found in time 
filters for signal analysis [20]. 

In more detail, motivated by the recent realisation of an atom laser in a waveguide 
[21], we shall focus on the longitudinal beam features and explore the space-time 
evolution of effective one dimensional sources represented by "source" boundary 
conditions of the form 

i;{x^O,t)^x{t)e-'^'^\ Vi>0. (1) 

The aperture function modulates the wave amplitude at the source and is therefore 
responsible for the apodization in time, that is, the suppresion of the fringes in the beam 
profile. However we shall show that, for long enough pulses or for multiple overlapping 
pulses, such effect holds only in a given time domain because of the revival of the 
diffraction in time. 

Localised sources [22, 23] represented by an amplitude iIj{x = 0,t) ~ e"*'^"* have 
been used for understanding tunnelling dynamics or front propagation and related time 
quantities such as the tunnelling times [24] . Moreover, this approach has shed new light 
on transient effects in neutron optics [25, 26, 8] and diffraction of atoms both in time 
and space domains [27, 28]. The connection and equivalence between "source" boundary 
conditions, in which the state is specified at a single point at all times and the usual 
initial value problems, in which the state is specified for all points at a single time, was 
studied in [29]. Notice that we shall use the source conditions (1) for simplicity, but the 
pulses could also be examined with the initial value approach. 

In section 1 we shall review the sudden switching on of the source since any other 
case is a combination of this elementary dynamics. Single pulse formation with standard 
apodizing functions is examined in sec. 2, and we describe the dynamics for a smoothly 
switched-on source in sec. 3; Section 4 is devoted to the multiple pulse case treated with 
different aperture functions, and the paper ends with a discussion. 

1. Turning-on the source 

We shall consider the dynamics of a source of frequency ujq = hkQ/2m, which is turned on 
in a horizontal waveguide. This has the advantage of cancelling the effect of gravity off. 
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avoiding the decrement of the de Broghe wavelength X^b — 2Tr/ko due to the downward 
acceleration [21]. Note that the effects of an external linear potential on the diffraction 
in time phenomenon can also be taken into account following [30]. At the same time 
the dynamics is effectively one-dimensional provided that ujq -C uj± where uj± is the 
transverse frequency of the waveguide. 

The sudden approximation for turning on a source is well-known [1, 23, 29, 24] 
but briefly reviewed here for completeness. The aperture function is then taken to be 
the Heaviside step function, Xo(^) = ^{t)- The inverse Fourier transform of the source 
iIjo{x = 0, ko, t) = e~^°*e{t) is given by 

ilj{x = 0,ko,uj) = F-^['iIjo]{x = 0,ko,uj) = / i^oi^ = 0, ko,t)e'''*dt 

(2) 



^/2^^ u} -uiQ + iO 
For x,t > 0, the wave function evolves according to 



i /"°° 
ipo{x,ko,t) = — duj 
27r y_oo 



. . oo ikx—iiot 
6 / - C- 



u — ujo + iO^ 

where 1mA; > 0. This can be written in the complex fc-plane by deforming the contour 
of integration to r_|_ which goes from — oo to oo passing above the poles. 



il:o{x,ko,t) = 7^ / dk2k 
27r Jr^ 



1.2 _ h.2 



fdk (^ + ^)e^'-'^. 



271" Jr+ \k + ko k — ko^ 
Since one of the definitions of the Moshinsky function M{x, k, t) is precisely 

■u' _-k''^t 

M(.,k,t) = -l^,k'^^, (3) 

(The Moshinsky function can be related to the complementary error function, see 
Appendix A.), the resulting state can then be simply written as 

V^o {x, ko, t) = M (x, ko, —)+M (x, -/Co, — ) . (4) 
\ m J \ m J 

Such combination of Moshinsky functions is ubiquitous when working with quantum 
sources and will appear throughout the paper. The source density profile corresponding 
to Eq.(4) exhibits the characteristic oscillations of diffraction in time [1, 2]. We 

are interested in describing more general switching procedures. However, given their 
dependence on pulse formation results, its discussion will be postponed to section 3. 



2. Single pulse 



In this section we study the formation of a single-pulse of duration r from a quantum 
source modulated according to a given aperture function x^^\^)-i the superscript denoting 
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the creation of a single pulse. Notice that some apodizing functions were already 
explored in the quantum domain [28], but we shall present a more general description. 

In particular, let us consider the family of aperture functions Xn\t) = 
sin"'(i7t)0(t)©(r — t) with Q = vr/r, see Fig. 1(a). The inverse Fourier transform 
of the source amplitude at the origin is given by 

1 2~"e**^'^~'^°)'^/^rr('l 



V'i')(x = 0,A;o, 



n] 



UJ 



27r. 



{2+n)TT+{u}—uio)T 

2-K 



{2+n)n—{ui—wo)T 
2n 



(5) 



where T{z) is the Gamma function 



The energy distribution of the associated 



wavefunction [29] is proportional to a;^/^|'0ri 



(1), 



X 



0,a;)p. As shown in Fig. 1(b), the 



smoother the aperture function the wider is the energy distribution, which will affect 
the spacetime profile of the resulting pulse. We shall next illustrate the dynamics for 
the n — 0, 1,2 cases. The solution for an arbitrary aperture function is presented in 
Appendix B. 

2.1. Rectangular aperture function 

For the case Xo^'*(^) — Q(t)Q{r — t) corresponding to a single-slit in time the dynamics 
of the wave function is given by 

V'?^ (x, ko, t- r) = V'o {x, ko,t)-e{t- T)e-^'^°^Vo {x, ko,t-r). (6) 

Note that if t < r, before the pulse has been fully formed, the problem reduces to that 
discussed in the previous section, being x{t) — ©(*)■ Equation (B.4) is also vahd for 
these times, in contrast to previous works restricted to t > r [25, 26, 27, 28]. The same 
will be true for the following pulses. 



2.2. Sine aperture function 

Let us consider now Xi^^i^) — sm(Qt) Q(t) ©(r — t). Then, defining uj± 
the corresponding wavenumbers k± — ^y2muJ±/h, 

^i\x, ko,t;T) = ^-'^ailj^^\x, k^^t-r). 



ujq^Q. and 
(7) 



2.3. Sine- square aperture function 

The sine-square (Banning) aperture function is given by 



x!^^\t)=sm\nt)Q{t)Q{r-t) = - 



1 - - cos{2Vtt) 



Q{t)Q{T-t). (8) 



Introducing kg = A/2m(c<;o ± 2Vi)/h with /3 = ±, leads to 



i^2\x,kQ,t]T) = - 



^o\x,ko,t]r) - ^^ip^^\x,k,3,t;r) 



/3=± 



(9) 
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Figure 1. a) Family Xn{t) of single- pulse aperture functions. The smoothness at 
the edges increases with in this order: rectangular (n = 0, continuous line), sine 
{n — 1, dashed line), sine-square (n = 2, dot-dashed line), b) The aperture function 
in the frequency domain, F~^[xn ^]{(^ — uio), broadens as a result of apodization, with 
increasing n. 



where the effect of the apodization is to substract to the pulse with the source momentum 
two other matter- wave trains associated with k/3, all of them formed with rectangular 
aperture functions. 

The effect of a short r on the energy distribution has already been studied in [2, 28] 
from which a time-energy uncertainty relation was inferred. Indeed, the uncertainty 
product reaches an approximate minimum of 2h if the pulse duration is shorter than 
To = tt/cjo- The space-time dynamics for a source apodized with a rectangular aperture 
function is illustrated in Fig. 2, where it is shown that for r > Tq the quantum average 
of the position operator follows the classical free-particle trajectory whereas if r < Tq 
the quantum pulse is sped up. Fig. 3 shows for different Xn{t) = 0, 1,2) the effect 
on the velocity distribution which is only centered at hko/m for r > Tq. Concerning 
the apodization, the higher the value of n the narrower is the mean width of Xn{t) and 
the larger the shift on the mean velocity. The smoothing of the Xn{t) is reflected on the 
suppression of the sidelobes in the velocity distribution. 

3. Smooth switching: apodization and diffraction in time 

Instead of forming a finite pulse as discussed in the previous section, some sources are 
prepared to reach a stationary regime with constant flux after the initial transient. The 
extreme case is the sudden, zero-time switching of Sec. 1. In this section we shall study 
how the matter-wave dynamics is modified for slow aperture functions with a given 
switching time r. In order to do so, we introduce the family of switching functions 

{0, t <0 
Xit), < t < r 
1, t>T 

where it is assumed that x(t) = 1, and which is represented in Fig. (4)a for 
x{t) = sin^ flgt with n = 0,1,2 and Qs = vr/2r. Note that for n = 0, one recovers 
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Figure 2. (color online) Spacetime effects of the time-energy uncertainty relation. 
A pulse of duration r < Tq is sped up with respect to the source frequency. The 
spacetime dynamics is illustrated for a ^^Rb source with r = 1ms with hko/m = Icm/s 
(left column, Tq — 0.0467ms) and hko/m = O.lcm/s (right column, Tq = 4.67ms), 
modulated with different apodizing functions to form a single pulse. The dashed line 
reproduces the classical trajectory. The grey scale changes from dark to light as the 
function values increase. Figures in the same column show the suppression of the 
sidelobes in the probability density with increasing n. 



the sudden aperture Xo(^) = 0(^) which maximises the diffraction in time fringes. For 
n = 1 (n=2) the source amphtude follows half a sine (square)-lobe. The general result 
for the time-dependent wavefunction simply reads, 

^^(x, ko, t- r) = e-'"°"V^o (x, ko, t - r) e{t - r) + ^« (x, A^o, t; r), (10) 

where the superscript s stands for switching and tp^^ refers to a pulse of length r and 
apodizing function x- 

Figure 4(b) shows the effect of a finite switching time on the oscillatory pattern at a 
given time. As a result of the smooth aperture functions new frequencies are singled out, 
leading to a progressive suppression of diffraction in time. Moreover, for a given Xn{t) 
function, a larger switching time r increases the apodization of the source, washing out 
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Figure 3. (Color on-line) Distortion of the velocity distribution for pulse duration 
smaller than Tq. The corresponding pulses to Fig. 2 are plotted in velocity space 
for a ^''Rb source with r = 1ms apodized with different aperture functions and a) 
hko/m = Icm/s (Tq = 0.0467ms) and b) O.lcm/s (Tq — 4.67ms). The triangle marks 
the classical velocity. 
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Figure 4. a) Family Xn{t) of switching-on aperture functions. The smoothness at 
the edges increases with n. b) Apodization of the ^^Rb beam profiles with increasing 
smoothness of the aperture function (higher n), and fixed switching time Ts — 0.5ms 
at t = 20ms, and hk^)/m — 1 cm/s. 



the fringes in the probabihty density, see Fig. 5(a). (In this sense, the effect of a finite 
band- width source is tantamount to a time dependent modulation, as was shown in [23].) 
However, it is remarkable that the effect of the apodization is limited in time because of 
a revival of the diffraction in time. The intuitive explanation is that the intensity of the 
signal from the apodization "cap", see Fig. 6, decays with time whereas the intensity 
of the main signal (coming from the step excitation) remains constant. For sufficiently 
large times, the main signal, carrying its diffraction-in-time phenomenon, overwhelms 
the effect of the small cap. One can estimate the revival time by considering that the 
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Figure 5. Revival of tlie diffraction in time, a) The oscillations in the beam 
profile characteristic of the diffraction in time are gradually suppressed with increasing 
switching time r. The beam profile is shown at t = 20ms, and hko/m — Icm/s. For 
increasing r the amplitude of the oscillations diminishes and the signal is delayed, b) 
During the time evolution there is a revival of the diffraction in time which is again 
clearly visible for times greater than tr — ujqt^ = 1ms, hko/m = 0.5cm/s, and 
tr — 16.8ms). Both a) and b) refer to a ^^Rb source which is switched on following a 
sine-square function. 
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Figure 6. Aperture function decomposition. In a switching proccess the aperture 
function has a finite transient component responsible for the apodization (shaded area) . 
Such apodization cap is followed by a step excitation for t > t which leads eventually 
to the revival of diffraction in time. 



momentum width of the apodization cap will be of the order 

Ap ^ ^2mh{{jjQ -f Q^) — \/2vrLh{ijjQ — Vt^) ~ ^/2mhuJo- , (11) 

and taking into account the, essentially linear with time , spread of the cap half-pulse 
described by tp^^ in Eq. (10) (we have checked numerically that the classical dispersion 
relation Ax{t) ~ Aa;(0) -|- Aptr/m holds). 

The ratio of the spatial widths corresponding to the initial and evolved cap ip^^ is 
given by 

^ _ Ax(0) Ax(0) 

Ax{tr) Ax{0) + Aptr/m' ^ ' 
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Taking Ax{Q) ~ ^/2hJ^JmT and imposing R — 1/2 leads to the expression for a 
revival time scale 

tr ^ ouor^, (13) 

which is analogous to the Rayleigh distance of classical diffraction theory but in the 
time domain [31]. The smoothing effect of the apodizing cap cannot hold for times 
much longer than tr- This is shown in Fig. 5(b), where the initially apodized beam 
profile develops eventually spatial fringes, which reach the maximum value associated 
with the sudden switching at t ~ 2tr- 



4. Multiple pulses 

We show next that the single pulse aperture functions x^^K^): described in section 2 
are useful to consider matter wave sources periodically modulated in time. Indeed, the 

A^-pulse aperture function can be written as a convolution of X^^K't) with a grating 
function. In particular, if the grating function g{t) is chosen to be a finite Dirac comb 
gN{t) = '^j'=o 5{t — jT), it is readily found that 

„ N-l 

X^''\t)=x^'\t)*g^{t) = / dt'x^'\t-t')g^{t') = Y^ x^^)(t-jT), (14) 

j=o 

which describes the formation of N consecutive x^^^-pulses, each of them of duration r 
and with an "emission rate" (number of pulses per unit time) 1/T. 

If we are to describe an atom beam periodically modulated with such apodizing 
function, the space-time evolution simply reads, 

iV-l 

)(a:, ko, t) = J2 ^o, t - jT; T)e{t - jT) (15) 

j=0 

being a linear combination of the single-pulse wavefunction ^*^^)(a;,t) studied in section 
2. The Q{t — jT) corresponds to the condition that the j-th pulse starts to emerge only 
after jT. 

A similar mathematical framework can be employed to describe an atom laser in 
the no- interaction limit. The experimental prescription for an atom laser depends on the 
outcoupling mechanism. Here we shall focus on the scheme described in [12, 13]. The 
pulses are obtained from a condensate trapped in a MOT, and a well-defined momentum 
is imparted on each of them at the instant of their creation through a stimulated Raman 
process. The result is that each of the pulses tpx^ does not have a memory phase, the 
wavefunction describing the coherent atom laser being then 

N-l 

<l>f\x,ko,t) = ^4^)(x,A;o,f-jT;T)e(i-iT), (16) 

j=0 

compare with (15). 

It is interesting to impose a relation between the outcoupling period T and the 
kinetic energy imparted to each pulse hjo, such that ujqT = Itt, where if / is chosen an 



Time m,odulation of atom sources 



10 



0.03 
0.02 
0.01 
0.00 

0.004 
0.002 

0.000 

0.015 
0.010 
0.005 
0.000 







a) 








iiii 


iiiiiiv^ 


b) 

/V 








c) 



50 



100 



150 



0.03 
0.02 
0.01 
0.00 

8e-05 
4e-05 
Oe+00 
0.003 
0.002 
0.001 
0.000 
0.06 
0.04 
0.02 
0.00 



x(|j,m) 



d) 



e): 



V 



200 



400 



600 



x(|j,m) 



Figure 7. Beam density profile of a ^'^Na atom laser with r = 0.833ms, T = O.Sms, 
in the coherently constructive ujqT = IOtt (5.9 cm/s) (a), destructive ujqT = IItt 
(5.9 cm/s) (b), and incoherent case (c), with same conditions than a) with a sine- 
squarc aperture function. The conditions equivalent to the experiment reported in [12] 
(r = 0.83ms, T = 0.05ms) are studied for the constructive (d), destructive (e) and 
incoherent case (f). If the pulse is not apodized the beam profile becomes noisy as 
shown in (g) for rectangular %, in the same conditions than (d). In all cases the norm 
is relative to the incoherent case. 



even (odd) integer the interference in the beam profile ■^^^^ is constructive (destructive) , 
see Fig. 7. However, in the case of a periodically chopped atom beam ipx'\ Eq. (15), 
the additional phase e'"*'^"-'^ singles out the constructive interference independently of 
the parity of 

So far we have assumed that the source is coherent. The interference pattern 
between adjacent pulses is affected by the presence of noise in the phases [32]. The 
average over many realisations in which the phase aj of each pulse varies randomly in 
[0, 27r] leads to the density profile. 



N-l 

e^"^^^\x, ko, t - jT; r)0(t - jT) 

j=0 



2v 



(17) 



which reproduces the incoherent sum of the pulses, p{x,t) = '^j\ilJn\x,ko,t — 

jT-T)\'e{t-jT). 

One advantage of employing stimulated Raman pulses is that any desired fraction 
of atoms can be extracted from the condensate. Let us define s as the number of atoms 
outcoupled in a single pulse, ipx^ ■ The total norm for the A'^-pulse incoherent atom laser 
is simply N'inc = sN. A remarkable fact for coherent sources is that the total number 
of atoms A^c outcoupled from the Bose-Einstein condensate reservoir depends on the 
nature of the interference [33]. Therefore, in Figs 7 and 8 we shall consider the signal 
relative to the incoherent case, namely, / Afinc- 

Note the difference in the time evolution between coherently constructive and 
incoherent beams in Fig. 8: At short times both present a desirable saturation in the 
probability density of the beam. However, in the former case a revival of the diffraction 
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Figure 8. Time evolution of a ^■^Na atom laser beam modulated by a sine-square 
function with r = 0.833ms, T = 0.05ms, moving at 5.9cm/s in the a) coherently 
constructive, b) incoherent case, and c) rectangular single pulse approximation to a) 
(the snapshot a.t t ~ 200ms is omitted for clarity). Note that Xd = hkot/m is the 
classical trajectory. In all cases the norm is relative to the incoherent case, see (17) 
and the paragraph bellow. 
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Figure 9. (color online) Many-pulse aperture function X2^''(^) with Teff = 5t and 
increasing number of pulses N: 5, (solid line), 10 (dashed line), and 20 (dot-dashed 
line). In the highly overlapping regime X2^\t) can be approximated by a a single 
pulse rectangular aperture function of duration r^ff — {N — l)T + t. 



in time phenomena takes place in a time scale ujot"^, similar to the smooth switching 
considered in section 3, whereas the later develops a bell-shape profile. 

Interestingly, whenever the overlap amongst different pulses is coherently 
constructive or with phase memory (as in the periodic chopping of a beam), an effective 
single-pulse aperture function can be introduced. Indeed, if r/T S> 1 as in the actual 
experiments [12, 13], it follows that Xeffit) = J2f=i xi^ ~ jT) oc 9(t)0(re// — t) for any 
X, where the effective pulse duration is Tg// = (A^ — 1)T -|- r. More precisely it reaches a 
plateau in the time interval [r, nT — t] but has both a front and a back cap of duration 
r associated with the switching on and off. This is graphically illustrated in Fig. 9 
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for X2^\'t)i describing the formation of N pulses with a sine-square modulation X2^^(0- 
Moreover, the quality of this approximation improves as the time of evolution becomes 
greater than the revival time, see Fig 8(c). Clearly this is only valid in the highly- 
overlapping regime whereas if the overlap is restricted to few pulses no saturation of the 
beam profile is observed. This is an important simplification to optimise the aperture 
function and design the beam profile of an atom laser. 



5. Discussion 



We have studied the dynamics of non-interacting boson sources which are modulated in 
time with arbitrary aperture functions. Starting with different switching-on procedures 
we have accounted first for single pulse formation and systematically described the 
multiple pulse case. Whenever the modulation is sudden, the spreading of the source 
exhibits oscillation in the density profile, which is the taletelling sign of the quantum 
diffraction in time phenomenon. For apodized sources with smooth aperture functions 
the oscillations are suppressed, thus bringing to the quantum domain the classical 
apodization techniques usual in signal analysis and Fourier optics. However, for long 
or multiple overlapping pulses the smoothing is limited by a time scale tr, associated 
with a revival of the diffraction-in-time phenomenon. Multiple pulse formation has been 
examined for different possible cases. For constructive or for random averaged phases, a 
saturation effect occurs providing a fiat intensity profile before t,., which later develops 
spatial fringes and a bell-shape profile respectively. For both types of beam formation 
the initial fiat beam region is smooth for apodized pulses but noisy for square ones. 
While these results are already quite relevant for an understanding of pulse formation 
and intended applications of the atom laser, the next natural step will be to extend the 
present analysis to systems with atom- atom interaction. 
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Appendix A. The Moshinsky function 

The Moshinsky function can be related to the Faddeyeva w{z) and complementary error 
function [34, 35], 

M(x,k,T) :^ ——w{-z), where 2; = — /r ( /c I, (A.l) 

2 2 V r/ 
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and w{z) is defined as 

_ 2 If e""^ 

w(z) := e exici-iz) = — du , (A. 2) 

ITT Jr_ u — z 

r_ being a contour in the complex 2;-plane which goes from — oo to oo passing below 
the pole. 

Appendix B. Arbitrciry aperture function 

In what follows we will be interested in pulses modulated with an arbitrary aperture 
function x(t) which is zero outside the interval [0,r]. An analytical result can be 
obtained in such case by means of Fourier series. Noting that 



^(^t)= y2 Cre'^e{t)e{T-t) with c, = - / x{t)e-'^dt,{B.l) 

r=-oo ^ 

the wavefunction at the origin can be conveniently written for all times as an infinite 
series 

oo oo 

i^^ix = 0, A;o, t)^ c^e-<'"°+^)*e(i)e(T - = ^ c^e-^'^'-*e(i)e(T - t), (B.2) 

r=— oo r=— oo 

which evolves for x, t > as 

oo 

^'x(^)^0,^;t) = Y Cri^''Q\x,kr,t]T) (B.3) 
r=— oo 

where 

'il)'i\x, kr, t; r) = Vo {x, K, t) - e{t - r)e~^'^'-^'0o {x, K, t - r) (B.4) 



and hkr — y^jnhJr-, with the branch cut taken along the negative imaginary axis of 
ujr — uJo-\- I'KrjT. This is an amusing expression for it relates the dynamics of a source 
with an arbitrary aperture function with that of rectangular pulses of different coherent 
sources. As a caveat, less general decompositions may be simpler as in the yti'it) case 
discussed in section 2. 
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